Vortex formation in a fast rotating Bose-Einstein condensate 
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We study rotational motion of an interacting atomic Bose-Einstein condensate confined in a 
quadratic-plus-quartic potential. We calculate the lowest energy surface mode frequency and show 
that a symmetric trapped (harmonic and quartic) Bose-Einstein condensate breaks the rotational 
symmetry of the Hamiltonian when rotational frequency is greater than one-half of the lowest energy 
surface mode frequency. We argue that the formation of a vortex is not possible in a non-interacting 
as well as in an attractive Bose-Einstein condensate confined in a harmonic trap due to the absence 
of the spontaneous shape deformation, but it can occur which leads to the vortex formation if we add 
an additional quartic potential. Moreover, the spontaneous shape deformation and consequently the 
formation of a vortex in an attractive system depends on the strengths of the two-body interaction 
and the quartic potential. 



PACS numbers: 03.75.Lm, 05.30. Jp 



I. INTRODUCTION 



After the creation and realization of the Bose-Einstein 
condensate (BEC) in an alkali atomic trapped Bose g as , a 
new field, namely strongly correlated atomic system, has 
been opened up [1] . It is always interesting to study the 
consequences of an external perturbation to the system. 
One way to perturb the system is by applying external 
rotation to the system. The rotation of a macroscopic 
quantum fluid exhibits interesting counter-intuitive phe- 
nomena. For example, an atomic interacting(repulsive) 
BEC confined in a rotating harmonic trap produces quan- 
tized vortices for a sufficiently large rotation [2-4]. The 
theoretical calculation for critical rotation frequency (f2 c ) 
based on purely thermodynamic arguments [5] is signif- 
icantly smaller than the observed value of f2 c . For suffi- 
ciently high stirring frequencies of the trap, the system 
deforms spontaneously [6] to lower the energy barrier 
which prevents the nucleation of vortices at smaller f2 
[7] and this spontaneous deformation is detected at the 
ENS experiment which leads to the vortex nucleation [8] . 
In fact, vortices are nucleated when the stirring potential 
with low anisotropics resonantly excites the quadrupole 
mode of the condensate, and induces large amplitudes 
oscillations of the condensate, resulting in a dynamical 
instability [9]. 

The harmonically trapped BEC becomes singular 
when the rotation frequency is equal to or greater than 
the harmonic trap frequency, because the outward cen- 
trifugal force counteracts the inward force from the har- 
monic trap. We can eliminate the singularity at O > 1 
if we consider an additional stiffer radial potential (say, 
quartic potential for simplicity). There is a growing in- 
terest about the effect of an anharmonic potential on 
the properties of a rotating BEC [10-14]. Recently, the 
quadratic-plus-quartic potential has been achieved exper- 
imentally at ENS by superimposing a blue detuned laser 



beam to the magnetic harmonic trap [15]. 

In spite of the fact that the Q c in a harmonically 
trapped BEC calculated from the thermodynamic ar- 
guments does not match with the experimental values, 
many authors have studied f2 c in a BEC confined in 
harmonic-plus-quartic trap based on the thermodynamic 
arguments [11-14]. In this work, we argue that a vor- 
tex formation is possible in an attractive as well as non- 
interacting Bose gas confined in a quadratic-plus-quartic 
potential and calculate the critical rotational frequencies 
based on the spontaneous shape deformation of rotating 
Bose gas. 

This paper is organised as follows. In section II, 
we write down the model Hamiltonian and variational 
macroscopic wave function, from which we calculate the 
variational energy. By minimizing the variational energy, 
we derive three coupled equations of the variational pa- 
rameters. In section III, we calculate the lowest energy 
surface mode frequency of an interacting BEC confined 
in a quadratic-plus-quartic potential. In section IV, we 
discuss how spontaneous shape deformation transfer an- 
gular momentum and vortex nucleation occurs in a har- 
monically trapped BEC. In section V, we discuss the pos- 
sibility of the shape deformation and vortex nucleation of 
a BEC trapped in a quadratic-plus-quartic potential. We 
present summary and conclusions of our work in section 
VI. 



II. THE SYSTEM AND THE HAMILTONIAN 

The equation of motion of the condensate wave func- 
tion is described by the mean-field Gross-Pitaevskii equa- 
tion, 



(i) 



i 



where V{f) = \m[(u 2 x 2 + uj 2 y 2 + ^ 2 z 2 ) + \uj 2 (x 2 + y 2 ) 2 ] 
is the quadratic-plus-quartic potential with the small an- 
harmonic term A. Here, we assume uu 2 = cu 2 (l + e) and 

uj 2 = loq(1 — e), where e = ^fq^f) is a very small defor- 
mation parameter due to the rotating deformed poten- 
4Z rr?~ ^ s tne strength of the mean-field 
xp y — yp x is the ^-component of the 



tial. Also, g = 
interaction, L z 
angular momentum operator, and Qq is the trap rota- 
tion frequency. When e = 0, it preserves the circular 
symmetry of the system. 

For simplicity, we shall study quasi-2D system. The 
wave function in the z-direction is separable and is given 

by, 



1>(z) = 



(0^) 1/2 



(2) 



where a 2 = After integrating out the z-component, 

we get the effective 2D Gross-Pitaevskii energy functional 
with the modified interaction strength g 2 = 2\ // 2nhuj z aa z 
[16]. One can write down the Lagrangian density corre- 
sponding to the quasi-2D system as follows: 



ih . , dip* 



dt' 



(3) 



+ (^IW>| 2 + V(x,yM 2 + |H 4 - QpL z ij>). 

Here, we use the time-dependent variational method 
[17] to study the properties of a rotating BEC confined 
in a quadratic-plus-quartic potential. In order to obtain 
the evolution of the condensate we assume the most gen- 
eral Gaussian wave function, 

1>(X,Y,t) = C( t ) e -iK*)* 2 +/3(*)y 2 -2i7(t)^] 5 (4) 



where C{t) = J N ^P is the normalization constant, and 



D = ol\0\ — 7j\ The new co-ordinates X and Y are 
the dimcnsionlcss variables, X = — , Y = — where 

a = is the oscillator length and 2lo 2 = uj 2 + lo 2 

is the mean frequency. Further, a(t) = oei(t) + ia 2 {t), 
0{t) = 0i{t) + i0 2 (t) and 7(f) = 71 (t) + ij 2 {t) arc the 
time-dependent dimcnsionlcss complex variational pa- 
rameters. The ai and 0i are inverse square of the con- 
densate widths in x and y directions, respectively. 

We obtain the variational Lagrangian L by substitut- 
ing Eq. (4) into Eq. (3) and integrating the Lagrangian 
density over the space co-ordinates, 

1 
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27172) + {on +0i)D 

+ (a 2 +7 2 )/3 1 + (/3 2 2 +7 2 2 )ai 

- 2(a 2 + /3 2 )7i 72 + ( ai + p\) + e(p\ - ai) 

A (3a 2 + 2aiA + 3/3 2 + 4 7l 2 ) 

2 D 
+ PD 3 ' 2 + 2Q{j 1 {0 2 - a 2 ) + 72 ( ai - /3i))], (5) 



where P = 2 A /- ( - jv 1 - )a is the dimcnsionlcss parameter 

that indicates the strength of the mean-field interaction 
and it can be positive or negative depending on the sign 
of the s-wave scattering length a [18]. 

The variational energy of the rotating condensate at 
equilibrium is given in terms of the inverse square width 
of the condensate along the x and y directions and the 
phase parameter <5 72 , 



^ = i[(a 1 +/? 1 )+7 2 2 (^- 



1 . . 1 1 
Pi a\ pi 



1 1 A 3 2 
+ e{ ~ 1 ' 0\ ) + 2X + «iA 

+ p^p\ + 2n l2 (±--— )], 

Pi ai 
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2) 



(6) 



where f2 = Oo/^o- We have used the fact that a 2 = f3 2 = 
71 = 0, and also 72 = 5j 2 to obtain the above equation. 

One can get the equilibrium value of the variational 
parameters, a±o, 0io and 720 by minimizing the energy 
with respect to the variational parameters, 




a( (l 



/3 2 o(l + T7 



and 



) = 7 2 - 2^720 + 1 + e ■ 



W 3 1 ^ 
«10 P10 

(7) 




7 2 2 + 2^720 + 1 



720 = 



{0io ~ aio) 
{P10 + «io) 



e + A( + — , 

"io P10 

(8) 



(9) 



The Eqs. (7), (8) and (9) describes how the shape of 
the condensate changes due to the rotation. The relation 
between the phase parameter and the variational widths 
is also obtained in [6,19]. The above equation can be 
re- written as 



1 



< ar > 

< y 2 > 



0w 
«io 



(» + 720) 

(ft -720)' 



(10) 



This implies that > 720 since 77 is always positive. 
The average angular momentum per particle is given 

by, 



< L z > = 1 J_ J_. 

Nh 2 720l /3 10 a w - 



» (1-r?) 2 
2 77(0:10 + /3io) ' 



(11) 



This relation explicitly shows how the angular momen- 
tum transfers to a trapped BEC due to the spontaneous 
shape deformation and consequently how the vortices ap- 
pears in the rotating BEC by calculating the average an- 
gular momentum per particle [20] . We would like to esti- 
mate the critical angular frequency from Eq.(ll) by using 
the constraint < L z >= Nh. 
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III. LOWEST ENERGY SURFACE MODE 
FREQUENCY 

To calculate the quadrupole mode frequency of BEC 
confined in a quadratic-plus-quartic potential, we expand 
the Lagrangian in the following way: a = a w + 5a\, 
= 0io + 50i > an d 7 = <$7i • We keep only the second 
order deviations from their equilibrium values. Then the 
Lagrangian quadratic in the deviations takes the follow- 
ing form: 



_ 1 5a 1 5a 2 5(3x502 . ^7i^72 
— 4^^2 1 ~ r z ~ 
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(12) 



For simplicity, first we set il = 0. Using the Euler- 
Lagrangian equation of motion, we can get the following 
two coupled equations of 5a\ and S0\. 
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(13) 



(14) 
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For an isotropic trap, e = 0, and we set Sai = —50\ to 
calculate the quadrupole mode (m z ± 2) frequency which 
is given by, 



where Rq is the equilibrium radius of the system, and this 
can be obtained from the real solution of the cubic equa- 
tion: (1 + -j)i?Q — i? ^ 4A = 0. The above quadrupole 
frequency u q is valid for all interaction (repulsive, at- 
tractive) strength with a small anharmonic term (A) . For 
A = 0, the critical value of the interaction strength (above 
which the system is unstable) is P = —2. For A ^ 0, and 
P < —2, there is no real solution and the system col- 
lapses. Note that the quartic potential does not affect 
the critical strength above which the system is unstable 
compared to the harmonic trapped BEC. The quadrupole 
mode frequency uj q vs. P is shown in the Fig.l. 
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FIG. 1. Quadrupole mode frequency tu q vs. interaction pa- 
rameter P for A = 0.005 (solid) and A = 0.05 (dashed-dot). 

For A = = P, it reproduces the known result 
for the quadrupole mode frequency of a harmonically 
trapped non-interacting BEC. Also, for A = 0, it gives 
the known result of the quadrupole mode frequency of a 
harmonically trapped interacting BEC [21]. This result 
shows that the Lu q is decreased compared to the non- 
interacting case due to the effect of the repulsive inter- 
action. Similarly, the Lu q is also increased compared the 
non-interacting case due to the effect of the attractive in- 
teraction. But small anisotropy (A ^ 0) increases the uj q 
drastically in all cases. Note that, the u q depends on two 
parameters, A and P. For a harmonically trapped BEC, 
uj q always decreases when P increases, whereas for the 
additional quartic term, u q first decreases for increasing 
P, but after some value of P, it starts increasing with P. 
The value of P at which u q starts increasing depends on 
A. For example, P c = 40 for A = 0.005 and P c = 10 for 
A = 0.05. 

The centrifugal term — QL Z shifts the quadrupole mode 
frequency by ±2f2. The lowest energy surface mode fre- 
quency is given by, 



^ q = (4+ 2 -^-PRl), 



(15) 



UJ-2 



4 + f -P* 



2n. 



(16) 
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IV. ROTATING BEC CONFINED IN A 
HARMONIC TRAP 

Using the three coupled equations (7), (8) and (9) for 
the three variational parameters, we get the following 
equation for the ratio of the widths: 

(1 + e - ft 2 )?/ 4 + 2(1 + e - ft 2 )?/ 3 + 2e?/ 2 (17) 
+ 2(ft 2 - 1 + e)ry + (ft 2 - 1 + e) + ^V^K 1 + e - fi2 )^ 3 
+ (1 + 3e - 5ft 2 )?? 2 + (5ft 2 + 3e - 1> + (ft 2 — 1 + e)] = 0. 

The above relation for the aspect ratio of the system is 
valid for all interaction strength including repulsive as 
well as attractive system also. The same kind of relation 
for the phase parameter 72 is obtained in [19]. 

First, we consider a symmetric case (e = 0) and P > 0. 
One would expect naively that the ratio of the square of 
the widths, 77 becomes one for all rotation frequency. But 
this is not the case for all rotation frequencies. 7/ is al- 
ways one when < ft < u> q /2, where u> q is the quadrupole 
mode frequency which is given by 

u q = sjl-PRl (18) 

But r\ has three solutions when ft is greater than ^f. 
One solution is still one and the other two solutions are 
less than one and greater than one. But the system has 
the lowest energy when 77 is either less than or greater 
than one. Now we discuss how the angular momentum 
is transferred to the system by the quadrupolar defor- 
mation. When the rotation frequency is in the interval 
of < ft < u) q /2, 77 is always one and it is clear from 
Eq. (11) that the average angular momentum per parti- 
cle is zero, although we are rotating the system. In other 
words, the system does not responds to the external rota- 
tion till ft = ^ . This is a signature of the irrotationality 
of the condensate. Once the system starts deforming, 
it starts transferring angular momentum (see Eq. 11) 
and vortex starts nucleate when the average angular mo- 
mentum per particle is one. The existence of the critical 
angular frequency at which the angular momentum per 
particle suddenly goes up from to 1 is the direct mani- 
festation of the supcrfluid of the condensate. The average 
angular momentum of each particle has been measured 
experimentally which is one when a single vortex is nu- 
cleated [3]. One can do the same analysis for e ^ 0. For 
large P and very small e, u q ~ \/2, then a single vortex 
will form when ft ~ 1/V2 <~ 0.7 which is experimentally 
observed. 

For non-interacting system, set P = in Eq. (17), 
then one can get the following quartic equation for 77: 

(1 + e - ft 2 )?? 4 + 2(1 + e - ft 2 )?? 3 + 2e?/ 2 
+ 2(ft 2 - 1 + e)r]+ (ft 2 - 1 + e) = 0. (19) 



In the non- interacting system, 77 is complex if the rotat- 
ing frequency ft is in the interval [ yl — e, yl + e] . Only 
one real solution exists for ft < \f\ — e and ft > y/1 + e 
[6,19]. For symmetric trap potential (e = 0) it is clear 
that there is only one real solution (77 = 1) exist for any 
rotating frequency ft. The lowest energy surface mode 
frequency is 2loq and hence instability will occur when 
ft = Do at which the system is destabilized. There is no 
spontaneous shape deformation and hence one can not 
produce a vortex in the non-interacting system because 
the angular momentum per particle is always zero. We 
have also checked that for finite but small e, 77 changes 
very slowly with ft which is not sufficient to transfer one 
unit of angular momentum per particle to the system and 
hence there is no vortex formation in the non-interacting 
system even for finite e. 

Now we consider an attractive BEC. We will always 
consider |P| < 2, since the condensate will be in stable 
state if \P\ < 2. In this attractive case, the quadrupole 
mode frequency is given below. 



^ = V 4+ 2^rk- (20) 

For any values of \P\ < 2, uj q > 2. It clearly shows and 
confirmed by the solution of the Eq. (17) that the shape 
deformation will not take place in the range of ft = to 
ft = 1. In other words, 77 is always one. Since there is no 
shape deformation and consequently vortex nucleation is 
not possible. Note that authors in Ref. [22,23] showed 
that the vortex formation in a rotating attractive BEC 
in harmonic trap is prohibited due to the fact that all 
the angular momentum is carried by the center of mass 
motion. Our analysis also shows why there is no vortex 
formation in an attractive BEC in harmonic trap. One 
can do the same analysis for e 7^ 0, but the deforma- 
tion is not enough to transfer unit quantum of angular 
momentum to each particle. 

V. ROTATING BEC CONFINED IN A 
QUADRATIC-PLUS-QUARTIC POTENTIAL 

In this section we will discuss the effect of the quar- 
tic potential on the possibility of a vortex formation in 
a harmonically trapped BEC. Using the three coupled 
equations (7), (8) and (9), we get the following two cou- 
pled polynomial equations: 

= (a 3 0/ 3 3 + 2a 4 /? 2 + a 5 w (3 w )(l + ^J^) 

2 V Qio 

- a? /3io(3ft 2 + 1 + e) + (a 10 (3 3 w + 2a 2 /? 2 )(ft 2 - 1 - e) 

- A(5a 2 /3io + 7aio/3 2 + a 10 + 3/3 3 ) (21) 

and 
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= (a? /3? + 2a? # + a 10 # )(l + f y^g) 

- a 10 /3? (3« 2 + 1 + e) + (a 3 w /3 10 + 2a? /3 1 2 )(O 2 - 1 - e) 

- \(5a 10 f3 2 10 + 7a 2 w /3 10 + 3a? + /3 3 w ). (22) 

To show that there is a circular symmetry breaking of a 
BEC confined in a quadratic-plus-quartic potential, we 
put e = and consider large P (say P = 1000) and solve 
those two coupled equations. We find that there is a cir- 
cular symmetry breaking when the rotating frequency is 
greater than one-half of the quadrupole mode frequency 
u) q . We plot r] vs. SI for e = and A = 0.005 which is 
given in Fig. 2. In this case, = 0.887179 and r\ is always 
one for fi = to Q = 0.887179. When Q > 0.887179, r) 
has three real and positive solutions and r\ < 1 or rj > 1 
corresponds to the lowest energy state. As the system 
starts deforming itself, it transfers angular momentum 
to the system. When fi = 0.905, each particle gets an 
unit angular momentum, i.e. < L z >= Nh. The vortex 
formation can occur for SI < 1 as well as SI > 1, de- 
pending on the parameter P and A. For A = 0.05, the 
critical rotation frequencies are Sl c = 1.2, 1.422, 1.554 for 
P = 100, 500, 1000, respectively. 
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FIG. 2. The variation of r\ vs. SI for A = 0.005 and large P. 

The shape deformation also occurs for non-interacting 
and attractive Bose gases. To discuss about the shape de- 
formation in non-interacting system, first we set P = 
in the equations (21) and (22). We find that the rj has 
three real and positive solutions when Q > like in 
the condensate in the harmonic trap. For quantitative 
discussion, we take A = 0.005, then the system starts 
deformation when SI > = 1.0123. In the previous sec- 
tion we have seen that the spontaneous deformation takes 
place due to the two-body interaction. But here spon- 
taneous shape deformation exists in an ideal Bose sys- 
tem confined in a quadratic-plus-quartic potential. The 
average angular momentum per particle is always zero 



when SI = to SI = 1.0123, but < L z >= Nh at 
SI = 1.027. We give another example for A = 0.05. For 
this anharmonic strength, the system will starts deform- 
ing when SI > ^j- = 1.10895 and the vortex will form 
when < L z >= Nh at SI = 1.21. It shows that as we 
increase the strength of the anharmonic term, the criti- 
cal rotational frequency is also increasing. Note that the 
vortex formation occurs only when SI > 1 because ^ 
is always greater than 1. The variation of r\ vs. SI for 
non-interacting BEC is shown in the Fig. 3. 
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FIG. 3. The variation of r? vs. Q, for P = 0, A = 0.005 and 
e = 0. 

It is worth to point out that the shape deformation 
can even occur in the non-interacting Bose gas if it is 
confined in a quadratic-plus-quartic potential. 

In the previous section, we have seen that harmonically 
trapped attractive BEC can not hold a vortex. This is 
because of the absence of the shape deformation in the 
rotating system. Now we will see what happens if we 
add the quartic potential in Hamiltonian. Once we add 
the quartic potential, we can access a large rotation fre- 
quency (SI > 1) for which the rotating attractive sys- 
tem is not destabilized. We assume P = —0.015 and 
A = 0.01 and then we get %^ = 1-02606 for the sym- 
metric trap. By solving the equations (21) and (22), we 
see that the spontaneous shape deformation starts when 
SI <~ 1.02606 and transferring an angular momentum to 
the system. Each particle gets an unit angular momen- 
tum, < L z >= Nh, when SI = 1.052. We give another 
example. Assume, P = —1.5 and A = 0.05 and then we 
get — 1.66215 for the symmetric trap. By solving the 
equations (21) and (22), we get three real and positive 
solution of -q when SI > . But the existence of the three 
solutions of rj does not mean that the system will starts 
deforming itself. Actually, the lowest energy state can 
be obtained only for r/ = 1. It means there is no shape 
deformation for these parameters and hence a vortex can 
not appear under this conditions. In this attractive sys- 
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tcm, the shape deformation can occur which leads to the 
vortex nucleation, depending on the strengths of the two- 
body interaction and the quartic potential. When there 
is no shape deformation for a given large P and A, the 
center of mass motion is excited [22,23]. These results 
are also consistent with the recent results of Lundh et 
al. [14]. Note that ^ is always greater than 1, so that 
deformation takes place only when > 1. 



VI. SUMMARY AND CONCLUSIONS 



In this paper we have calculated the quadrupole mode 
frequency. We have shown that when the rotating fre- 
quency of an interacting BEC confined in the quadratic- 
plus-quartic potential is greater than one-half of the 
quadrupole mode frequency, the system itself starts de- 
forming from a circular shape to an elliptic shape and 
hence it breaks the rotational symmetry of the Hamilto- 
nian. We have shown how this self-deformation transfer 
an angular momentum to the rotating system and con- 
sequently creates a vortex. 

We have argued that the formation of the vortex is 
not possible due to the absence of the spontaneous shape 
deformation of a harmonically trapped rotating attrac- 
tive BEC as well as an ideal BEC. But the shape defor- 
mation and consequently the vortex formation can oc- 
cur in a rotating attractive BEC if it is confined by an 
additional quartic potential. We have also shown that 
the shape deformation occurs not only due to the two- 
body repulsive interaction, it can occur even in the non- 
interacting as well as attractive BEC if it is confined 
in a quadratic-plus-quartic potential and consequently it 
produces a vortex. However, the shape deformation and 
consequently forming a quantized vortex in a rotating at- 
tractive BEC depends on the strengths of the two-body 
interaction and the quartic potential. 

Note that, in actual experiments the lowest energy 
surface mode is excited through the rotating potential 
with small anisotropics. In this work we have studied 
(for simplicity) the vortex formation in a symmetrically 
trapped BEC, but we can easily extend our work with 
small anisotropy e. 
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